High-precision determination of the electric and magnetic form factors of the proton 
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New precise results of a measurement of tiie elastic electron-proton scattering cross section per- 
formed at the Mainz Microtron MAMI are presented. About 1400 cross sections were measured 
with negative four-momentum transfers squared up to Q^ = 1 (GeV/c)^ with statistical errors be- 
low 0.2%. The electric and magnetic form factors of the proton were extracted by fits of a large 
variety of form factor models directly to the cross sections. The form factors show some features 
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at the scale of the pion cloud. The charge and magnetic radii are determined to be {te)^ = 

0.879(5)stat.(4)syst.(2)model(4)groupfm aud {rli)^ = 0.777(13)stat.(9)syst.(5)inodel(2)groupfm. 

PACS numbers: 13.40.Gp, 14.20.Dh , 25.30.Bf 
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The proton is composed of relativistic light con- 
stituents, i.e. quarks and gluons, as manifest in "deeply 
inelastic" lepton scattering experiments. On the other 
hand, its global charge and magnetization distributions 
are contained in the form factors of elastic electron scat- 
tering. These provide a very instructive connection be- 
tween models based on effective degrees of freedom and 
the QCD-based parton theory of the nucleon [Ij] . A par- 
ticularly intriguing question is the existence of a direct 
signal for a pion cloud in the form factor hypothesized on 
the basis of pre-2003 data in Ref. [3]. Since Yukawa the 
virtual pion is considered as the mediator of the nucleon- 
nucleon force. Today it is viewed as an effective degree 
of freedom of the nucleon originating from the sponta- 
neous breaking of chiral symmetry producing the God- 
stone bosons of QCD. An exciting question, addressed in 
this letter, is whether it could be seen directly. 

A quantity particularly sensitive to the existence of a 
pion cloud is the proton charge radius. The determina- 
tions from electronic hydrogen Lamb shift measurements 
[3, 3] agree with those from e-p scattering if Coulomb 
and relativistic corrections are applied [5|-l7|. This work 
presents a novel determination of the electric and mag- 
netic form factors and radii using a direct fit of form 
factor models to the electron scattering cross sections. 

The Mainz accelerator MAMI provides a cw electron 
beam with energies up to 1600 Me V with very small halo 
and good energy definition. However, in this experiment 
only electron energies up to 855 MeV were used since 
the higher energies were not yet available at the time 
of the experiment. Taking advantage of the three high- 



resolution spectrometers of the Al collaboration it was 
possible to measure the elastic electron-proton scattering 
cross section with a statistical precision of better than 
0.2% and extract the form factors up to a negative four- 
momentum transfer squared of Q^ — 0.6 (GeV/c)^. 

About 1400 cross sections were measured at beam en- 
ergies of 180, 315, 450, 585, 720, and 855 MeV covering 
Q-^ from 0.004 to 1 (GeV/c)^. In order to achieve high 
accuracy the experiment aimed to maximize the redun- 
dancy of the data. For the angular scans the spectro- 
meter angles were varied only in small steps so that the 
same scattering angle is measured up to 4 times with dif- 
ferent regions of the spectrometer acceptance, and parts 
of the angular range were measured with two spectrome- 
ters. The broad range of beam currents from below 1 nA 
to more than 10 /xA, a consequence of the large range 
of cross sections, required special attention to the deter- 
mination of the luminosity. Therefore, the current was 
measured redundantly with a fluxgate magnetometer and 
with a pA-meter connected to a collimator just down- 
stream the electron source. Furthermore, the relative lu- 
minosity was measured at all times with one of the three 
spectrometers at a fixed scattering angle. 

In lowest order, the elastic electron-proton scattering 
cross section is described by the Rosenbluth formula 
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where Ge and Gm sue the electric and magnetic Sachs 
form factors, rup is the proton mass, r = Q'^ /{Am'ic'^) and 

e = (l + 2 (1 + r) tan2(e'/2)) "^ with the electron scatter- 



ing angle 9. However, also electromagnetic processes of 
higher order contribute to the measured cross section, 
such as multiple photon exchange, vacuum polarization, 
vertex corrections, and the radiation of a real photon 
from the electron (Bethe-Heitler) or the proton (Born). 
The code simulating the cross-section integration over 
the acceptance includes these processes following the de- 
scription of ref. [8| which gives results compatible with 
ref . [9J . Our approach extends this by an explicit calcu- 
lation of the Feynman graphs of the Bethe-Heitler and 
Born processes on the event level. The simulation uses 
the standard dipole parametrization 
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as a sufficient approximation for the true form factors 
(/ip is the proton's magnetic moment divided by the nu- 
clear magneton). The division of the measured num- 
ber of elastically scattered electrons by luminosity and 
simulated acceptance-integrated (standard dipole) cross- 
section yields the measured normalized cross section; this 
procedure accounts for the radiative processes. 

Furthermore, Coulomb corrections according to 
ref. [101 have been applied, but no correction for the ex- 
change of two hard photons (TPE) in the scattering pro- 
cess since no unique prescription exists yet. 

The energy of the elastically scattered electrons Eout is 
shifted by the recoil energy of the proton. The finite reso- 
lution of the spectrometers, the external energy-loss pro- 
cesses, and the internal bremsstrahlung widen the peak 
in the AE = i5out,caic. - -^out.meas. spcctrum. The elastic 
events are selected by a cut around this peak. 

The targets used were 2 and 5 cm long cells filled with 
liquid hydrogen with walls made of 10 /zm thick havar 
foil. The primary source of background is the elastic 
and quasielastic scattering off the nuclei in these walls. 
As verified with empty-cell measurements, inelastic peaks 
are either small or outside the cut region. The amplitude 
of its simulated shape was fitted to the measured AE 
spectrum together with the simulated hydrogen peak. As 
shown in Fig.[T] the simulated shapes describe the mea- 
sured spectra very well. Therefore, after subtraction of 
the background from the raw spectra, the result is very 
insensitive to a variation of the cut around the elastic 
peak. In particular it shows the validity of the radiative 
corrections. The background contribution reaches up to 
10% but is below 4% for most of the data. 

For fixed Q^ Eq. ([T]) can be rewritten as an equation of 
a straight line in e in which G'^ and G^ are fit parame- 
ters. This "Rosenbluth" method is model independent to 
first order in the photon propagator. However, one has 
to take the data at constant Q^ for a sufficiently broad 
range of e. For a given energy and angular range this lim- 
its severely the kinematical range covered. This restric- 
tion can be avoided by a direct least-squares fit of models 
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FIG. 1. Typical AE spectrum. Red histogram: measured 
spectra (height of elastic peak: 5500), blue line: simulated 
background, green shaded band: data minus simulated hy- 
drogen peak and simulated background. The width of the 
band denotes the statistical uncertainty of the data. Data 
were measured at Eq = 450 MeV, Q^ = 0.045(GeV/c)^ 



for Ge and Gm to the measured cross sections for all Q^ 
and 6. Of course, it is mandatory to select a wide range 
of different models and test for any model dependence. 
In this work we have used single-dipole and sum-of-two- 
dipoles models, the phenomenological parametrization 
of ref. [2j , several variations of polynomial, reciprocal of 
polynomial, and spline based models, and a variant of 
the extended Gari-Kriimpelmann model [ll'|. A detailed 
study of the model dependence and the analysis method 
was performed using pseudo data generated according to 
a variety of previous parametrizations [2, \l2, uM ■ This 
study showed that all models with enough flexibility, i.e. 
polynomials, their reciprocals, and splines, are able to 
reproduce the input with an error smaller than the other 
systematic errors. 

Since a determination of the absolute normalization 
of the measurement to better than 1% is not possi- 
ble, the normalizations of the individual cross-section 
data sets are left as free parameters with the constraint 
Ge{m){Q'^ — 0) = l(/^p)- The fitted normalizations are 
well within the estimated uncertainty of 4% and have 
almost no dependency on the model. For the flexible 
models the spread between the largest and smallest nor- 
malization factor is below 0.3%. 

The results of the direct flts are compatible with the 
results from a classic Rosenbluth separation where such a 
comparison is possible. However, we flnd that the Rosen- 
bluth approach is more sensitive to systematic deviations 
and is therefore a less robust estimator of G'^ and G^^. 
However, the Rosenbluth separation allows us to check 
for deviations from a straight line caused by possible 
problems in the data or by higher-order processes like 
TPE. At the level of the uncertainty of the measurements 
no systematic deviations from straight lines were found. 

As the direct flts of models are nonlinear, standard er- 
ror estimation techniques for the flt are not guaranteed 



to be exact. Therefore, the confidence bands were calcu- 
lated with the Monte Carlo technique including the er- 
rors of the normalizations. We find that Monte Carlo and 
the linearization used in standard error propagation yield 
almost identical results for all but one model. The con- 
fidence bands presented here are the widely used point- 
wise bands, meaning that one expects the true curve to 
be with 68% probability within the band at any given 
single Q^, but not necessarily at all Q^ simultaneously. 
The Monte Carlo approach also allows one to construct 
simultaneous bands meaning that with 68% probability 
the true curve does not leave the band for the full range 
of Q^. It is somewhat involved to treat this problem 
with standard analytical methods IJ] . The simultaneous 
bands can be obtained from the pointwise bands shown 
here by scaling the latter by a factor of around 2.3 for 
the Q2 range up to 0.6 (GeV/c)^. 

The form factors extracted with the flexible models 
agree among each other to better than 0.25% in the 
Q2 range up to 0.5 {GeV/cf. They all fit the data 
equally well with x^/d.o.f. « 1.14 for d.o.f. « 1400. 
However, including the less flexible models one obtains 
1.16 ^ x^/d.o.f. ^ 1.29 and the agreement is only bet- 
ter than 0.6%. In Fig. [2] the results of the spline model 
for Gb, Gm and their ratio are shown, together with 
previous measurements and fits. The error bars of the 
previous data shown for Ge and Gm are statistical only, 
normalization uncertainties are typically of the order of 
a few percent. Since TPE corrections are not applied to 
any of the data, the corresponding non-TPE-corrected fit 
of Ref. [IS] is shown. In the plot of the ratio the fit to 
the TPE-corrected data of ref. [13[ is also included. 

The results for Ge exhibit a large negative slope rel- 
ative to the standard dipole at Q^ w giving rise to 
the significantly larger charge radius. This slope lev- 
els out around 0.1 (GeV/c)^ and remains constant up to 
0.55 (GeV/c)^ when the slope again becomes larger. In 
that region, however, only measurements at large scatter- 
ing angles for only two beam energies contribute so that 
the fit becomes less reliable and more sensitive to system- 
atic errors such as the neglect of TPE. For even higher 
Q^ measurements have been taken only at one energy 
and a separation of Ge and Gm is not possible. In the 
close-up for Ge there is an indication of a bump around 
0.15 (GeV/c)^, however, at the limit of significance. 

The magnetic form factor Gm deviates from earlier 
measurements. This may be related to the normalization 
at Q^ — 7- ignoring the wiggle seen by this experiment. 
The maximum and the minimum of the wiggle structure 
depend, of course, on the parameter of the dipole form. 
Also, it is not clear whether the older experiments include 
the proton contribution to the radiative corrections. 

The structure at small Q^ seen in both form factors 
corresponds to the scale of the pion of about Q^ « mj « 
0.02 (GeV/c)^ and may be indicative of the influence of 
the pion cloud [Ij]. 
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FIG. 2. Tlie form factors Ge and Gm normalized to the stan- 
dard dipole and Ge/Gm as a function of Q . Blaclc line: best 
fit to the data, blue area: statistical 68% pointwise confidence 
band, light blue area: experimental systematic error, green 
outer band: variation of the Coulomb correction by ±50%. 
The difi'erent data points depict previous measurements, for 
Refs. see 0, [l3i "^^ added the data points of [Isl-KLTj. Dashed 
lines are previous fits to the old data in [2, [l3] . 



While the deviation of Gm from previous measure- 
ments seems surprising at first glance, it reconciles the 
form factor ratios from experiments with unpolarized 
electrons, like this one, with those found with polar- 
ized electrons, especially with the high-precision mea- 
surements in ref. [17| . The previous Ge and Gm data 
are basically not compatible with the polarized measure- 
ments even when TPE corrections are applied. New re- 
sults from Jefferson Laboratory [l8|, |l9| with uncertain- 
ties of about 2% confirm this statement and are in excel- 
lent agreement with this experiment. 

The charge and magnetic rms-radii are given by 
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In the study of the model dependency through simu- 
lated data only the flexible models reproduce the input 
radii reliably. In the fits to the measured data the models 
can be divided into two groups: Those based on splines 
with varying degree of the basis polynomial and num- 
ber of support points and those composed of polynomials 
with varying orders. For the charge radius the weighted 
averages of the two groups differ by 0.008 fm. 

For the spline group we obtain the values 

{rlY = 0.875(5)stat.(4)syst.(2)modclfm, 

{tIjY = 0.775(12),tat.(9)syst.(4)modcifni 
and for the polynomial group 

{rlY = 0.883(5),tat.(5)syst.(3)modoifm, 

{rl,Y = 0.778(ti5)«tat.(10)sy«t.(6)modcifm. 

Despite detailed studies the cause of the difference be- 
tween the two model groups could not be found. There- 
fore, we give as the final result the average of the two 
values with an additional uncertainty of half of the dif- 
ference: 



0.879(5)stat.(4)syst.(2)modol(4)group 



fm, 



i^ll)^ = 0.777(13)stat.(9)syst.(5)modol(2)groupfm. 

These radii have to be taken with the applied correc- 
tions in mind. While the Coulomb correction used is 
compatible with other studies [5|, |6| a more sophisticated 
theoretical calculation may affect the results slightly. 

The electric radius is in complete agreement with the 
CODATA06 [13 value of 0.8768(69) fm based mostly on 
atomic measurements. It is also in complete accord with 
the old Mainz result '2l| when the Coulomb corrections 
p, |6| are applied. However, the results from very recent 
Lamb shift measurements on muonic hydrogen 22| are 
0.04 fm smaller, i.e. 5 standard deviations. This differ- 
ence is unexplained yet. The calculation of the Lamb 



shift in muonic hydrogen requires the solution of a rela- 
tivistic bound state problem (see Ref. |23| and references 
therein) . The deviation may be due to the distorted wave 
functions, significantly more distorted than in electronic 
hydrogen, necessitating the consideration of multiphoton 
exchange. 

The magnetic radius has a larger error than the charge 
radius since the experiment is less sensitive to Gm at 
low Q^. Its value is smaller than results of previous fits, 
however, it is in good agreement with ref. 2J|, who found 
0.778(29) fm from hyperfine splitting in hydrogen. 

The consequences of the results presented here for our 
picture of the proton are discussed in refjl[ . A full ac- 
count of this work will be published 25|, |23| . 
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